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nuclei of 59/2- and /7/2 subshell are studied using a schematic model inter- 



o 
o 

^^ ' action with four force parameters proposed recently. It is shown that the 

0^ ' model interaction well reproduces observed physical quantities: the double 

differences of binding energies, symmetry energy, Wigner energy, odd-even 

^ ' mass difference and separation energy, which testifies to the reliability of the 



model interaction and its p— n interactions. First of all, the double differences 



p\ ' of binding energies are used for probing the p-n interactions. The analysis 



reveals different contributions of the isoscalar and isovector p — n pairing in- 
teractions to two types of double differences of binding energies, and also 
indicates the importance of a unique form of isoscalar p — n pairing force with 
all J components. Next, it is shown that this p-n force is closely related to 
the symmetry energy and the Wigner energy. Other calculations demonstrate 
the significant roles of p — n interactions in the odd-even mass difference and 
in the separation energy at N = Z. 

PACS number(s): 21.10.Dr;21.10.Hw;21.60.-n;21.60.Cs 



I. INTRODUCTION 

With the advent of radioactive nuclear beams, the properties of nuclei beyond the proton 
stability line have attracted experimental and theoretical attention in recent years. Special 
interest is devoted to a unique aspect originating from the fact that protons and neutrons 
occupy the same orbits in nuclei with N ^ Z{ see [Q for review ). Consequently, one 
expects a strong proton-neutron {p — n) interaction because of the large spatial overlaps 
between proton and neutron single-particle wave functions. The correlation energies related 
to the p — n interaction have been extracted from the experimental binding energies |0-§[ • 
A double difference of binding energies has been analyzed with the aim of providing input 
to semiempirical mass formulas |@-3, and with relation to the clustering of nucleons as 
elementary modes of excitation in nuclei [^. It has recently been used for study of the p-n 
interactions [§,0, and discussed in terms of schematic and realistic shell model calculations 
p. This approach using the double difference of binding energies may provide details of the 
p — n interactions, the isoscalar (r = 0) and isovector {t = 1) p — n interactions. The p — n 
part of the isovector pairing correlation near N = Z has been studied in terms of algebraic 



model P] and compared with shell model Monte Carlo calculations ||TO|,|TT|. On the other 
hand, there are many discussions of the roles of the isoscalar pairing interactions (see Refs. 
|I^-P^ for instance). 



The experimental data indicate [y-^ that the symmetry energy accompanied by the 
so-called Wigner energy behaves according to the T(T + 1) dependence (T = \Tz\). This 
form could come from the isospin-invariant Hamiltonian. It has been recently proposed 
that the Wigner energy originates in the isoscalar pairing interaction |l^JT3|] . They pointed 
out that the Wigner energy cannot be solely explained in terms of correlations between 
the J = 1 isoscalar p — n pairs, and the isoscalar p — n pairs with the other J contribute 
significantly |T^. In fact, a recent shell model calculation |T^ with the J = isovector and 



J = 1 isoscalar pairing forces in the N ^ Z fp shell nuclei cannot explain the magnitude of 



the Wigner energy. On the other hand, another paper |T^ has discussed the fact that the 



degeneracy of the T = and T = 1 states in odd-odd nuclei with N = Z is produced by a 
balance of the symmetry energy and the J = isovector pairing correlation. The lowering of 
the T = states in N = Z odd-odd nuclei, according to our investigation p0|j2l[] , is caused 



by a unique form of isoscalar {p — n) pairing force including all J components. This result 
is consistent with that of Satula et al. [T^IT^ The uniqueisoscalar p — n pairing force, which 



can be expressed in a simple form including the T{T + 1) term, manifests a close relation to 
the symmetry energy. We shall discuss this matter by a concrete calculation in this paper. 
The odd-even mass difference (OEMD), the extra binding energy of a nucleus relative 
to its neighbors, is known to be an obvious experimental evidence of the pairing corre- 
lation [|22|. The pairing phenomena are well understood in terms of the proton-proton 
{p ~ p) or neutron-neutron {n — n) pair condensate, and described by the Bardeen-Cooper- 
Schriefer(BCS) theory p3|. The OEMD is often interpreted as a measure of the pairing 
gap (following the relation 12/A^/^ on the average) in medium-heavy and heavy nuclei. The 
OEMD displays, however a different feature in N ^ Z nuclei such as a special increase at 
N = Z. An ordinary estimation of the neutron or proton pairing gap from the OEMD is 
not applicable to these nuclei. On the other hand, it has recently been discussed that the 
OEMD in light nuclei is affected by deformation as well as J = pairing correlation P^-EB 



A further investigation of the OEMD should be made in N ^ Z nuclei. We shall discuss 
the influence of the p — n interactions on the OEMD. 

The development of recent radioactive nuclear beams facilities provides unstable nuclei 
beyond the line of proton stability. Experimental and theoretical investigations of proton 
emitters are increasing. Such phenomena allow a test of the various models on the proton- 
rich side. For N k, Z nuclei, one can expect that the p — n interactions also influence the 
separation energy. Since the p — n force is considered to be attractive, it might increase 
the separation energy. In fact, the calculated separation energies by all models without the 
p — n interaction are smaller than those of experiments at the N = Z nuclei. 

We need a reliable effective interaction to study the nuclear properties mentioned above. 
We have proposed an extension of the P + QQ model with four forces |^ , which reproduces 



quite well the experimental binding energies and energy spectra in N ^ Z nuclei of gQ/2 and 
fj/2 subshells. This model interaction including different types oi p — n forces is very suitable 
for our purpose to study various aspects of the p — n interactions. The main purpose of this 
paper is to study the p — n interactions, analyzing the double differences of binding energies, 
and to check the validity of our model examining various quantities such as the symmetry 
energy, the Wigner energy, the odd-even mass difference and the separation energy in nuclei 
near N = Z. 

The paper is organized as follows. In Sec. II, we first review our model proposed in the 
previous paper. Sec. Ill contains the analysis of the double differences of binding energies to 
probe the p — n interactions. The symmetry energy and the Wigner energy are discussed in 
Sec. IV. In Sec. V, the odd-even mass differences are analyzed in detail, and the two-proton 
separation energies are calculated in Sec. VI. Finally, Sec. VII gives conclusions. 

II. MODEL INTERACTION 

We have proposed the following effective interaction extended from the P + QQ force 



which is composed of four isospin-invariant forces (see Ref. [^ in detail): 



H = H,^ + Mnt (1) 

■"sp = 2^ ^aCapCap, [■^) 
ap 

Vint = v{Po) + v{QQ) + v{P2) + v;=', (3) 

viPj) = -\gjY.Y.PlMiA<^h)Y^PjMUcd), (4) 

Mk a<b c<d 

v{QQ) = -\x ■■T.T.QlM)T.Q2,Acd) ■, (5) 

fJ- abp cdp' 
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PjMini^^b) =pj{ab)AlM^^{ab), 
Qlf^piab) = q{ab)Bl^p{ab), 
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A'jMrMb)= Y. iJam^J^m,\JM)j:{7;P7;P'\r'^)^^, (9) 






Blpiah) = Y. {jam^Jbmp\2^i)ci^i-y^~"'^cpp, (10) 
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where po{ab) = J{2ja + l)Sab and p2{ab)=q{ab) = {a\\r'^Y2\\b) / y/E. We use the notation JM 
and TK, for the spin and isospin of a nucleon pair, respectively. The subscript p denotes the 
z components of isospin ±|. We also use the notation p = vr for a proton and p = u for a 
neutron. 

Here Hsp is a single-particle Hamiltonian and V^int contains the four forces: V{Po) stands 
for the isovector monopole pairing force, V{QQ) for the isoscalar quadrupole-quadrupole 
force, ^(-Pj) for the isovector quadrupole pairing force and VJJ"^ for the J-independent 
isoscalar p — n force. The first two forces in the interaction (0) are an extension of the 
conventional P + QQ force to the isospin-invariant one. The p-n part of the monopole and 
quadrupole pairing forces as well as the p — n component of the quadrupole-quadrupole force 
would be important for N k. Z nuclei. The last p — n force is very important for reproducing 
the experimental binding energy. It is important to note that V^j7° can be expressed as a 
simple form 

Kr = 4^°4(^ + l)-n, (11) 

where h denotes the total number operator of valence nucleons {h = hp + hn) and T is the 
total isospin operator. Our p — n interaction is composed of four different components and 
hence is useful in analyzing their respective contributions to various physical quantities. 
We applied the above hamiltonian to examine the binding energies and energy spectra 



of nuclei with A = 82 - 100 and A = 42 - 50 |2T|. We adopted only the 5(9/2 shell for 
nuclei with ^4=80-100 regarding the Z = N = 4Q core as inactive, and the fij2 shell for 
nuclei with ^4=40-50 regarding the Z = N = 20 core as inactive. It may be necessary to 
extend these model spaces for quantitative discussion. Our calculation itself indicates the 
insufficiency of the model space (/7/2)" about energy spectra. We used an extended model 
space (pi/2,5'9/2)" when comparing calculated energy levels with observed ones in nuclei 
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with A ^ 90, while experimental data near A = 80 seem to demand a further extension 
of the model space. The previous paper, however, has shown that the single j shell model 
is bearable for semiquantitative discussion about the nuclear binding energy. This simple 
model makes it possible to clearly see the roles of respective p — n interactions. We therefore 
employ the same single j shell model as that used in Ref. pT|, where the following force 
strengths are used: 

5-0 = 0.26, X = x{g(fl'9/2fl'9/2)}^ = 1-50, 

G2 = -g2{qi99/299/2)V = 0.35, k° = 0.925 in MeV for the c/9/2 shell region, (12) 

go = 0.59, X = x{q{f7/2f7/2)V = l.20, 

G2 = -g2{q{h/2h/2)Y = 0.90, A;° = 1.90 in MeV for the fr/2 shell region. (13) 

Our model with these sets of parameters is considered to be reliable for studying the p-n 
interactions in connection with physical quantities related to the binding energy in the (79/2 
and fi/2 shell nuclei. The mass dependence of the force parameter k^ is taken into account 
in some cases, but it does not change qualitatively the result. 

III. DOUBLE DIFFERENCES OF BINDING ENERGIES AND P-N 

INTERACTIONS 

We define the mth double difference of binding energies as follows: 

6V'^'^\Z,N) = 6^"'^B{Z,N), (14) 

where B{Z, N) is the nuclear binding energy. Here the operator 5^"^^ is defined as 

5('-)/(Z, N) = -—^[f{Z, N) - f{Z, N-m) 

-f{Z ~m,N) + f{Z -m,N -m)]. (15) 

The double difference of binding energies, SV^^\Z, N), was introduced for investigating the 
semiempirical mass formula &-ffl . This quantity is expected to roughly represent the p — n 



interactions between the last proton and neutron from the form of Eq. (pISl). Figure 1(a) 
shows the plot of 5V^^\Z, N) as a function oi A = N + Z for nuclei in the mass region y4=16- 



165. Experimental data are taken from Ref. ||2^. We see two separate groups in Fig. 1(a), 



namely, one is for the even- A nuclei (dots) and the other is for the odd-A nuclei (crosses). 
In both cases, shell effects at Z or A^ = 28,40,50,82 are present, while the patterns of 
dots and crosses are symmetric with respect to the average curve. It is now convenient to 
divide SV^^\Z, N) into two parts: the average part of the even-A and odd-A nuclei and the 
deviation from it. As seen in Fig. 1(a), the former is approximately written as Iq = 40/ A 
and the latter has opposite signs for the even-A and odd-A nuclei as follows: 

6V^'\Z,N)^Io + {-l)^h. (16) 



This expression was originally given by de-Shalit P,|28[| in the earliest investigations of the 



effective p-n interactions. Equation (|T6D describes the staggering with respect to the isotopes 
for even-y4 and odd-A nuclei. Large values of 6V^^\Z,N) for even-A nuclei (dots) near 
N = Z below A = 80 are notable. 

The data of 5V'''^\Z,N) are plotted in Fig. 1(b) as a function oi A = N + Z. (Our 
definition of ^V^^^-* has a sign opposite to that of Brenner's etal. 0|.) The values of SV^"^^ 
show a different behavior from 5V^^\ It is interesting that the staggering of SV^^'^ disappears 
in 5V^'^\Z, N). We see large scatters of dots and crosses for A < 80. These correspond to 
SV^'^'^ of nuclei in N ^ Z, and the values of SV^"^^ at N = Z are especially large. With 
decreasing mass A, SV^'^^ at N = Z increases. If one neglects the dots and crosses in N ^ Z 
nuclei, ^V^^^-* varies rather smoothly. This smooth trend is clear for A > 80 and continues up 
to heavy nuclei. This is due to the fact that there is no stable N k, Z nuclei with A > 80. 
Figure 1(b) clearly indicates the smooth systematic decrease of 5V'^'^^ with increasing mass 
A, which can be traced by the curve 40 /A. The deviations from the curve 40/A are small, 
and shell structure is not found. This general trend of 5V^'^'^ has long been known, and was 
discussed in several papers. In a recent paper [g], the dramatic spikes of 5V^'^'> at N = Z 
light nuclei were discussed in terms of both schematic and realistic shell model calculations. 



and the importance oi t = p — n interaction for the spikes was pointed out. The following 
relationship is derived from Eqs. ( p^ and ([T5|): 



+6V^^\Z -l,N) + 5V^^\Z - 1, iV - 1)]. (17) 

Substituting the empirical relationship ([T6|) into Eq. ([T7|) , in the large A limit we get 

5V'^^\Z,N)^lQ = A<d/A. (18) 

Thus, the systematic behavior of SV^'^\Z, N) for A > 80 can be explained from the relation 
(P^). ( Strictly speaking, there is a deviation from 40 /A due to the mass dependence of 
Iq.) Furthermore, it is obtained that 5V^"^\Z, N) for m = 3 — 6 have similar pattern to 
V^ \Z, N), and are also traced by the curve 40/^4. 

To analyze the double differences of binding energies, we now express the ground-state 
energy as follows: 

E{Z, N) =< H >= E,p + Ef»+««+^^ + E::' + E^^+,„ (19) 

^sp = (i^sp), (20) 

e::' = {v-'), (22) 

rpPo+QQ+P2 _ /-irPo+QQ+P2\ (r,o\ 

where ( ) denotes the expectation value with respect to the ground state. Here, \A^o+QQ+^2 
is the p — n parts of the Pq + QQ + P2 force, and V^^^^^^^ ^ is the proton-proton {p — p) 
and neutron- neutron {n — n) parts of the total interaction (Q). Figure 2 shows _E^o+QQ+^2 
and -E^^° as a function of valence-neutron number n„ for Nb and Mo isotopes. The p — n 
part of the PQ + QQ + P2 energy, £'^°+*^*5+'^2, exhibits a characteristic odd-even staggering in 
Nb isotopes, while -E^j7° gives a smooth line except for n„ = 1. On the other hand, for Mo 
isotopes ii;^o+QQ+-P2 varies smoothly as n„ increases, and indicates very different structure 
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from that of Nb isotopes. This can be attributed to extra energy for the odd-odd nuclei, 
which mainly comes from the t = 1 p — n part of the Pq + QQ + P2 force. 

Consider the double difference of ground-state energies, S^^^E{Z, N), using the operator 
5^^^ defined by Eq. (|15|). Since there are almost no contributions from the single particle 
energy Egp and Coulomb interaction to the double difference of binding energies as seen in 
the form of Eq. (|15|), one notices that 5^^^E(Z,N) is able to be compared directly with 
the experimental value SV^^\Z, N). Figure 3 shows the calculated and experimental double 
differences of binding energies, S^^'E{Z, N) and SV^^\Z, N), as a function of mass A = N+Z 
for the Nb, Mo, Tc, and Pd isotopes. The Nb and Mo nuclei near N = Z a.t the beginning 
of gg/2 shell region probably have the mixing of single particle levels, the pi/2, /5/2, and P3/2. 
As seen from Fig. 3(a)-(d), however, the agreement with experiments is quite good. Our 
calculation reproduces the staggering, and also predicts the highest spikes at N = Z nuclei 
though no experimental data are present. 



TABLE I. The components of 6^^'E{Z,N) for the Mo isotopes. The first and second columns 
denote the r = 1 components, and the third and fourth columns the r = components. 





T = 1 


r = 


Total 


N 


^(1) fpPo+QQ+P2 


S^''>E^^+,, 


SiDEQQ 5i^)El=^ 


5(i)i? 


42 


2.049 


0.001 


-0.001 -0.001 


2.050 


43 


-0.602 


-0.288 


0.451 0.464 


0.024 


44 


0.414 


0.398 


-0.232 0.460 


1.041 


45 


-0.511 


-0.096 


0.236 0.461 


0.089 


46 


0.346 


0.218 


-0.274 0.457 


0.747 


47 


-0.217 


-0.004 


-0.125 0.463 


0.117 


48 


0.018 


0.177 


-0.022 0.465 


0.638 


49 


-0.233 


0.000 


-0.186 0.461 


0.047 


50 


-0.014 


0.183 


-0.022 0.465 


0.611 



Let us now analyze the staggering and the highest spikes a,t N = Z . In Table I, the 



components of 6^^^E{Z,N) for the Mo isotopes are listed. The components S^^^ E^°'^'^'^~^^'^ , 
5^^^El^^ and 5^^^ E^t,^,^,^ are obtained using the definition (|15D of the operator 5^^^ for the 
respective parts of the ground state energies, Eqs. (^TD -(p3D. It is seen that the large value 
at N = Z = A2 comes from only the r = 1 component of §W E^l^'^^^'^^^ , and others are very 
small. Thus it is clear that the t = 1 p — n interaction of the Pq + QQ + P2 force is closely 
related to the large values of 6V^^^ at N = Z. On the other hand, 6^^^E{Z, N) foiN^Z 
exhibits staggering as seen in Table I [also see Fig. 3(b)]. The small values for odd A^ are 
due to the cancellation of r = 1 and r = components of 5^^^E(Z, N), and for even A^ both 
r = 1 and r = components contribute in phase. The value of S^^^E^^^ is for N = Z and 
^k^ for N > Z, though the tabulated values have numerical errors. 

Figure 4 shows 6'^'^^E{Z, N) and 6V'-'^\Z, N) as a function of mass A = N + Z hi the 
Mo, Tc, Pd, and Sn isotopes. The values of S^'^'E{Z,N) are a little bit smaller than the 
experimental ones but the agreement is quite well. Our calculation predicts large 5^'^^E(Z, N) 
at N = Z. The components of 5^'^^E(Z, N) are shown for the Mo isotopes in Table II. 



TABLE II. The components of 6'-'^^E{Z,N) for the Mo isotope. 





r = 


T = 1 


Total 


N 


Si2)EQQ 6^')El=' 


§m EP0+QQ+P2 


<5(2)i?,,+,. 


5^^)E 


42 


0.275 


0.694 


0.365 


-0.295 


1.039 


43 


0.237 


0.578 


0.128 


-0.125 


0.819 


44 


0.110 


0.462 


-0.094 


0.058 


0.536 


45 


-0.006 


0.462 


-0.083 


0.069 


0.441 


46 


-0.021 


0.461 


-0.083 


0.063 


0.419 


47 


-0.042 


0.461 


-0.092 


0.070 


0.397 


48 


-0.072 


0.463 


-0.099 


0.080 


0.373 


49 


-0.088 


0.463 


-0.110 


0.089 


0.355 


50 


-0.105 


0.464 


-0.123 


0.094 


0.331 



It is seen that 5^'^^E(Z, N) comes from only the t = p — n interaction, and the r = 1 com- 
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ponents are small because of the cancellation of S^'^^ E^°~^'^'^^^'^ {Z , N) and 5^'^'^Et,t^^i,„{Z, N). 
The r = component of QQ is small except for N = 42 and 43. (In other isotopes, this 
component is small for all A^, then this behavior for N = 42, 43 in the Mo isotopes is ex- 
ceptional.) The value of S^^^E^f is 3^74 for N=Z, 5A;78 for \N - Z\ = 1 and ky2 for 
\N — Z\ > 1, though the tabulated values have numerical errors. Thus SV^'^\Z, N) derived 
from the experimental binding energies are considered to be nearly attributed to the r = 
p — n force VJ^^^. Namely, we deduce an approximate relation 

6V^^\Z,N)^6^^^E{Z,N) 

^6^'^E::\Z,N). (24) 

This is consistent with the argument that SV^'^'{Z, N) vanishes if one neglects the r = 
p — n interaction in the shell model calculation with a surface 6 interaction for the 2s — Id 
shell §. 

Furthermore, we calculated 5^^^E{Z^N) and 5^'^^E{Z,N) for the I/7/2 shell nuclei using 
the shell model (/7/2)" with the force parameters (p!3D . Figure 5 shows the calculated and ex- 
perimental double differences of binding energies for the Ti and Cr isotopes. The agreement 
with experiments is very good. Our effective interaction (Pq + QQ + -P2 + K^;7°) reproduces 
well the experimental values of the double differences of binding energies also in the I/7/2 
shell region . This supports our model Hamiltonian being applicable to a wide range of 
nuclei. The good agreement tells us that the results in the gg/2 shell region are reliable and 
give good predictions. It should be also noted that the approximate relation Eq. (^) holds 
in this 1/7/2 shell region. 

It is now meaningful that the p — n correlation energy -E"^^" is expressed as 

Ei:' = -^k'i^i^ + 1) - TiT + 1)}, (25) 

for states with the total valence- nucleon number n and total isospin T from Eq. (^). For 
A^ — Z > 1, we can easily show that 

5(^)e;-° = y. (26) 

11 



The global behavior of SV^'^\Z, N) depending on 40/^4 as seen in Fig. 1(b) combined with 
the relations (p^ and (|26| ) suggests that in a wide-range view the force strength k^ might 
have 1/A dependence 

k' - f . (27) 

(Strictly speaking, since the double difference 5^'^^El'^^ with fc° = 80/^4 deviates from the 
curve 40 /y4, we need a higher-order correction with 1/A^^^ to reproduce the curve 40/A.) In 
fact, the force parameters k^ employed, 0.925 MeV for the lgQ/2 shell nuclei and 1.9 MeV 
for the 1/7/2 shell nuclei, reflect some A dependence. These values do not very deviate from 
the global value 80/A, if we compare them with the examples 80/A = 0.93 for A = 86 
and 80/A = 1.74 for A = 46. Certainly, if we impose the 1/A dependence on k^ like 
1.9 X (48/^4) in the calculations for the l/7/2-shell nuclei, the binding energies obtained for 
N ^ Z nuclei are reproduced better [^. This improves the double difference of binding 
energies 5E^'^\Z,N) as seen in Fig. 5, where the circles stand for the constant /c° and the 
crosses for the k'^ = 1.9 x (A8/A). Accordingly, the observed variation AO/A in 5V^'^\Z, N) 
is suggested to be mainly attributed to the global dependence 80/A on k^. The t = p — n 
force V^^^^ is possibly applicable to A^ > Z nuclei with T = Tz = {N — Z)/2. If it is true, 
the relation (|2^) holds for A^ > Z too. The p — n force VJJ"^ with the global parameter 
k^ = 80/A (with correction) may explain the smooth systematic behavior of SV^'^\Z, N) in 
the mass region A > 80 as seen in Fig. 1(b). This must be examined further. So far, we 
have not adopted the A dependence for the Pq, QQ, and P2 forces, because we do not have 
any strong demand to do so within the present calculations in a very tiny model space using 
a single j shell. It must, however, be necessary for our model when we make quantitative 
calculations in many j shells. 

IV. SYMMETRY ENERGY AND WIGNER ENERGY 

Let us next discuss the symmetry energy -Esym = cisjm{N — Z)'^ / A and Wigner energy Ew 
using the same set of parameters as the in previous sections. The experimental data indicate 

12 



that the symmetry energy accompanied by the Wigner energy is proportional to the T(T+1) 
where T = Tz = \N — Z\/2. Since -E^j7° includes the T{T + 1) term as seen in Eq. (p5D, 
both quantities must be closely related to the isoscalar p — n force V^J"^. Figure 6 shows the 
symmetry energy coefficient a{A) = 4asym in the expression Esyra+ Ew = a{A)T(T+l)/A for 
the /7/2 shell nuclei. The symmetry energy coefficient can be extracted by the treatment of 
Janecke and Comay |^^. We calculated the Coulomb-energy-corrected binding energies 



B* = B{exp) + _Ecoui(cal) following Caurier etal. [^. The calculated symmetry energy 
coefficients nicely reproduce the experimental data in Fig. 6. Where does the symmetry 
energy come from? We should now analyze the result obtained. If V{QQ) and V{P2) are 
eliminated from the total Hamiltonian, the Hamiltonian H^p + V{Pq) + V^^^ has SO (5) 
symmetry in the single-j shell approximation. The total energy of H^^ -|- V{Pq) -|- V^^^ in 
the single-j shell is specified by the total valence nucleon number n = rip -|- n„ and the total 
isospin T as follows 0: 

= en - -{^on(fi - ^) + fc°-(- + 1)} + -(^0 + Or(r + 1), (28) 

where E^^, Ep^ and E^^^ denote the expectation values of Hg^, V{Po) and V^j^^ with respect 
to ground states with n and T, respectively. From the coefficient of the T(T + 1) part in 
Eq. (p8D, the symmetry energy coefficient a{A) is expressed as 

~a{A) = ^{go + k')A, (29) 



which is proportional to the sum of the force strengths go and A;°. The parameter set (|T2[) 
gives the value d{A) = 1.2A5A in the /7/2 region. As shown in Fig. 6, the symmetry energy 
coefficient a{A) almost describes that obtained with the total Hamiltonian including V{QQ) 
and V{P2). Thus, we see that the symmetry energy in this region originates in the isoscalar 
p — n force V^^^ and isovector J = pairing force. Their contributions are 76% and 24%, 
respectively, in the present calculation. 

In the expression -Esym + Ew = a{A)T(T + 1)/A, the Wigner energy has the same 
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coefficient as the symmetry energy, and is expressed as Ew = a{A)T/A. Poves and Martinez- 
Pinedo pointed out that a shell model calculation with the J = isovector pairing force and 
J = 1 isoscalar pairing force in the N k, Z fp shell nuclei cannot explain the magnitude 



of the experimental Wigner energy [IS]. If we take the same parameter qq = 0.295 for 



the J = isovector pairing force as that of Ref. |T8|, the Wigner energy is estimated as 
Ew = go\N - Z\/4 = 3.54\N - Z\/A MeV for A = A8. This value is not very different 
from the result Ew = 3.04|A^ — Z\/A MeV they obtained. The empirical Wigner energy, 
Ew = 4:7\N-Z\/A MeV ^ or 37\N-Z\/A MeV [|1,|1] is very large compared with these 
values. Figure 6 tells that the Wigner energy cannot be reproduced without V^^^. If we 
introduce VJj"^, the Wigner energy becomes Ew = 37.4|A^ — Z\/A MeV for A = 60 from 
Eq. (p9D , which is consistent with the empirical formula. We can conclude that the isoscalar 
p — n force and isovector J = pairing force are origin of both the symmetry energy and 
Wigner energy. In particular, it is important to note that the isoscalar p — n interaction 
with all J components, not only J = 1, is crucial for reproducing the symmetry energy and 
the Wigner energy. 

V. ODD-EVEN MASS DIFFERENCE 

The odd-even mass difference in three-point and four-point expressions, 

A,{Z, N) = ^-^[B{Z, N+1)- 2B{Z, N) + B{Z, N - 1)], (30) 

A4(Z, N) = ^—^[B{Z, N + 1)- 3B{Z, N) + 3B{Z, N - 1) 

-B{Z,N-2)], (31) 

is often used to estimate the empirical pairing gap (for neutron) and to determine the 
pairing force strength. Figure 7 shows the experimental values of A3{Z,N) and A4{Z,N) 
as a function of A^ — Z in even-Z isotopes with proton number Z=20-30 and 36. In Fig. 
7(a), A^{Z, N) exhibits staggering around l.SMeV and has a notable peak ai N = Z, while 
A^lZ, N) has a hill near N = Z and N = Z + 1 but varies smoothly for N > Z + 2. The 
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OEMD, A3{Z,N) and A4{Z,N), are about 1.5 MeV on the average ior N > Z + 2. This 
value is usually regarded as a measure of the empirical neutron pairing gap. In addition, 
we notice the asymmetry of A3{Z,N) with respect to N — Z = 0. This may be due to 
the so-called Nolen-Schiffer anomaly [Q, an energy difference between neighboring mirror 
nuclei, which cannot be explained by the electromagnetic interaction. 

The calculated values of A3{Z,N) and A4{Z,N) are obtained by replacing B{Z,N) by 
the ground-state energy E{Z,N) in Eqs. (|30|) and (|3lD , since the Coulomb energy hardly 
contributes to A3{Z,N) and A4{Z,N). In our single j shell model, the single-particle 
energy has no contribution to A3{Z,N) and A4{Z,N). Figure 8 shows the calculated and 
experimental values of A^lZ, N) as a function oi N — Z for the Ca, Ti, Cr, and Fe isotopes. 
The agreement with experiments is quite good. Especially, the observed peaks at N = Z 
are well reproduced. 



TABLE III. The components of A3{Z,N) for Cr the isotopes. The first and second columns 
denote the r = 1 components, and the third and fourth columns the r = components. 





r = 1 


r = 


Total 


N 


AP0+QQ+P2 

'-^TTU 


AP0+QQ+P2 

'-^TrTT + UU 


aQQ 




A3(Z,Ar) 


21 


0.001 


1.553 


0.002 


0.000 


1.556 


22 


-0.027 


1.984 


0.106 


-0.004 


2.059 


23 


0.008 


1.657 


-0.107 


-0.003 


1.555 


24 


0.750 


1.500 


-0.009 


0.469 


2.710 


25 


0.003 


1.662 


-0.111 


-0.013 


1.541 


26 


-0.034 


2.007 


0.089 


-0.015 


2.047 


27 


0.001 


1.579 


-0.015 


-0.005 


1.560 



Let us analyze what interactions are important in A3(Z, A^). We separately calculated the 
contributions of the interaction energies E^l-^^^^^\ El=^ and E^lX^^^^^ to Ao,{Z, N), and 
denote them by A^o+^Q+^^ A;=° and Af°+^^+■^^ shortly. In Table III, the components of 
A3(Z, A^) are listed for the Cr isotopes. This table indicates the dominance of A^Jl^+iJi/^ ^) 
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namely the dominance of the hke-nucleon pairing correlations. The other components of 
A^lZjN) are very small except that the isovector parts of A^1~^^^^^^ and A^^° have large 
values at A^ = Z. We see their additional contributions for the large peaks of A3{Z,N) at 
A^ = Z in Fig. 8. The isovector and isoscalar p — n interactions are most cooperative with 
the p — p and n — n interactions in the N = Z nuclei. 

This situation is explained by illustrating the behavior of the respective interaction en- 
ergies in Fig. 9. The staggering of A3{Z,N) in Fig. 7 is almost attributed to that of 



E. 



P0+QQ+P2 

■KTT + UV 



in Fig. 9. The straight lines of the interaction energies i^^o+<9<9+-P2 and i?^^ ° go 



down as A^ increases and turn to the different directions at A^ = Z. The coincident "bends" 
at A^ = Z cause the increase (the peak) of A^{Z^N) according to the form of Eq. (pOD. 
These bends of the p — n interaction energies produce the increase of 5V^'^'> around N = Z. 
The bends give a special energy gain to the N = Z even-even nuclei, ^^Ti, '^^Cr, ^^Fe etc. 
The a-like four-nucleon correlations in these N = Z nuclei can be interpreted in terms of 
the characteristic behavior of the p — n interactions in cooperation with the like-nucleon 
interactions P^]. 

TABLE IV. The components of A^^+^^+^2 for the Cr isotopes. 



N 


aPo 

'-^TTTT+l^l' 


aQQ 

'-^TTTT + UU 


'-^wn+uu 


A-P0+QQ+P2 


41 


1.182 


0.374 


-0.003 


1.553 


42 


1.919 


0.589 


-0.524 


1.984 


43 


2.120 


0.279 


-0.742 


1.657 


44 


1.901 


0.398 


-0.799 


1.500 


45 


2.129 


0.281 


-0.748 


1.662 


46 


1.950 


0.594 


-0.537 


2.007 


47 


1.211 


0.381 


-0.013 


1.579 



According to Refs. [^,0 , on the other hand, the OEMD in light nuclei is strongly affected 
by deformation originated in the Jahn- Teller mechanism [RSI. It is interesting to see what 
interactions contribute to A^°^j;J^^ ^ being the main part of A3. Table IV presents the 
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respective contributions of the Pq, QQ and P2 forces to A^lXuif~^ ^ in the Cr isotopes. The 
dominant component is A^°_,_^^ as expected, i.e., about 2.0 MeV for | N — Z |> 2 and about 
1.2 MeV for N = Z ±3. In addition, there are considerably large contributions of ^^jf+uu 
and A^^_,_^^. Since the QQ correlation is intimately related to the nuclear deformation, 
the positive contribution of A^^^^ is consistent with the conclusion of Ref. |^. The 
contribution of quadrupole pairing force P2, however, is negative, and is larger than that 
of the QQ correlation for | N — Z \> 2. This is easily understood by the fact that the 
quadrupole pairing correlation breaks the J = Cooper pairs of neutrons. The present 
calculation tells us that the quadrupole pairing correlation probably cancels the effect of the 
QQ correlation or the deformation on the OEMD value p5|j26[| . 

Figure 10 shows calculated and experimental values of A^lZ, N) as a function of N — Z 
for the Ca, Ti, Cr, and Fe isotopes. The agreement with experiments is quite good except 
for the Ti isotopes. The calculation reproduces the hill near N = Z and Z + 1 and also 
the gentle behavior of A4{Z,N) near the value 1.6 MeV in the region of A^ > Z + 1 and 
N < Z. The increase of A4(Z, A^) at N = Z and N = Z + 1 is explained in terms of the 
same mechanism as that of A3{Z,N) at N = Z, which is caused by the coincident bends 
at A^ = Z of the two graphs in Fig. 9 that illustrate the variations of the p — n interaction 
energies. The isovector p-n interactions of the Pq + QQ + -P2 force and the isoscalar p — n 
force V^;7° are important at N = Z also for A4(Z, A^). 

VI. TWO-PROTON SEPARATION ENERGY 

We have investigated several quantities related to the nuclear binding energy in the 
previous sections. The calculated data include binding energies of nuclei close to the proton 
drip line. It is interesting to look at the two-proton separation energy, experimental data 
of which has been accumulated by the radioactive beam. It provides the possibility for 
studying new decay modes such as diproton emission. Some nuclei around ''^Ni are expected 
to possibly be two-proton emitters. There is a large deviation between theory and experiment 
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for the two-proton separation energy up to now. All the predictions by the Hartree-Fock 
(Bogoliubov) and relativistic Hartree-Fock (Bogoliubov) treatments p6| . |37[] underestimate 
the two-proton separation energies ai N = Z nuclei. This discrepancy could be due to 
the lack of p-n interactions in these treatments. As seen in the OEMD in Sec. V, the p-n 
correlations cooperate with the p — p and n — n correlations especially at N = Z nuclei. 
We can expect that the p — n interactions have a considerable influence on the two-proton 
separation energy. 

We calculated the two-proton separation energies S2p for the /7/2 shell nuclei with Z = 
20 — 28 and A^ = 20 — 28. In the calculation, the force strength kP is chosen as k^ = 
1.9 X [48 /A) which was used in the previous paper ||21]], because the 1/A dependence of k^ 
improves the binding energy and is also supported by the discussion in Sec. III. Calculated 
values of 5*2^ are compared with experimental data taken from Ref. |2^ in Fig. 11. Here we 



subtracted the Coulomb energy following Caurier etal. ||31[|. The agreement is good. The 
observed values of S2p are reproduced well ai N = Z nuclei. Again, the isoscalar p — n force 
V^J"° plays an important role in the two-proton separation energy at A^ = Z as well as in 
the other quantities discussed above. The model space (/7/2)"' and the set of parameters 
used are not appropriate for nuclei with large Z and A^, strictly speaking. According to the 
experience in this paper, however, the results on the different quantities of binding energies 
might be still meaningful. We note the prediction of our calculation that ^^Ni, ^^Co and 

2p 



^^Ni are possibly unstable, ^°Ni is bound, and S2p of "^^Fe and ^'^Co are close to zero. 



VII. CONCLUSION 

We have studied the p — n interactions using the functional effective interaction with four 
force parameters which reproduces the energy levels and binding energies oi N ^ Z nuclei 
considerably well. 

First, we analyzed the double differences of binding energies SV^^"^ and 5V^'^\ because 
the two quantities are expected to directly represent the p — n interactions. Our effective 



interaction reproduces fairly well the experimental values of SV^^^ and SV^'^\ and their 
characteristic behaviors, in the gg/2 and f-j/2 shell nuclei. The staggering of 6V^^^ is due to 
the competition between r = 1 and r = components. The large spike of 5V^^^ at N = Z 
is attributed to the r = I p — n interactions of the Pq + QQ + P2 foce, and that of 5V'^'^^ 
at N = Z contrarily represents the r = p — n interactions. The observed values of 5V'^'^^ 
with respect to mass A are approximated by the curve 40/ A. This curve may be explained 
by granting an A dependence on the r = p-n force strength /c°. 

Second, our effective interaction has also reproduced well the symmetry energy and the 
Wigner energy for the /7/2 shell nuclei. The strong t = p — n force V^j^^ with assistance of 
the r = 1, J = pairing force is important to explain the magnitudes of the two quantities 
in this region. It should be noted that the isospin parts of the two forces are proportional 
to T{T + 1) and their sum is directly related to the symmetry energy in the mass. The A 
dependence of the symmetry energy coefficient seems to be determined mainly by that of 

Third, our effective interaction has described well the observed values of the odd-even 
mass difference (A3 and A4) for the fi/2 shell nuclei. The cooperation of the p—n interactions 
with the like-nucleon interactions is remarkable at N = Z. It causes the rise of A3 and A4 
at N = Z. The characteristic behaviors of the p — n interaction energies at N = Z (see Fig. 
9) have an important effect not only on the double differences of binding energies but in the 
odd-even mass difference. 

We have briefly touched on the two-proton separation energy 6*2^ using the calculated 
binding energies. The calculation indicates a considerably large effect of the p — n force on 
5*2^ at N = Z. We noted the prediction of our calculation for 5*2^ near the /7/2 proton drip 
line. 

The present investigations have shown that the p — n interactions cause the notable 
behaviors of the observed quantities related to the binding energy near N = Z in nuclei 
where valence protons and neutrons occupy the same shells. Furthermore, our calculations 
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suggest that the p — n interactions (especially V^j^^) are important for describing these 
quantities over a wide range oi N > Z nuclei including the neutron drip line. 

All the results support the usefulness of the functional effective interaction composed of 
the four forces Pq, QQ, P2 and V^J^°, and clarify the essential roles of VJJ"^ in the physical 
quantities related to the binding energy. The present calculations, however, have been 
carried out in the single j shell model. Calculations in more realistic model spaces are in 
progress. 

Acknowledgments 

The authors are grateful to J.-Y. Zhang for helpful discussions. 



20 



REFERENCES 

[1] A. L. Goodman, Adv. Nucl. Phys. 11, 263(1979). 

[2] J. Janecke and H. Behrens, Phys. Rev. C 9, 1276(1974). 

[3] A. S. Jensen, P. G. Hansen, and B. Jonson, Nucl. Phys. A 431, 393(1984). 

[4] J. Janecke, E. Comay, NucL Phys. A 436, 108(1985). 

[5] G. G. Dussel, R. J. Liotta, and R. P. J. Perazzo, Nucl. Phys. A388, 606(1982). 

[6] J. -Y. Zhang, R. F. Casten, and D. S. Brenner, Phys. Lett. B 227, 1(1989). 

[7] G. Zaochun and Y. S. Chen, Phys. Rev. C 59, 735(1999). 

[8] D. S. Brenner, C. Wesselborg, R. F. Casten, D. D. Warner, and J. -Y. Zhang, Phys. 
Lett. B 243, 1(1990). 

[9] K. T. Hecht, Phys. Rev. B 139 (1965) 794; Nucl. Phys. A 102, 11(1967). 

[10] J. Engel, K. Langanke, and P. Vogel, Phys. Lett. B 389, 211(1996). 

[11] K. Langanke, D. J. Dean, S. E. Koonin, and P. B. Radha, Nucl. Phys. A 613, 253(1997). 

[12] W. Satula and R. Wyss, Phys. Lett. B 393, 1(1997). 

[13] W. Satula, D. J. Dean, J. Gary, S. Mizutori, and W. Nazarewicz, Phys. Lett. B 407, 
103(1997). 

[14] J. Engel, S. Pittel, M. Stoitsov, P. Vogel, and J. Dukelsky, Phys. Rev. C 55, 1781(1997). 

[15] D. Rudolph, et al., Phys. Rev. Lett. 76, 376(1996). 

[16] D. J. Dean, S. E. Koonin, K. Langanke, and P. B. Radha, Phys. Lett. B 399, 1(1997). 

[17] K. Kaneko and Jing-ye Zhang, Phys. Rev. C 57, 1732(1998). 

[18] A. Poves and G. Martmez-Pinedo, Phys. Lett. B 430, 203(1998). 



21 



[19] P. Vogel, |nucl-th/98050T5 . 



[20] K. Kaneko, M. Hasegawa, and Jing-ye Zhang, Phys. Rev. C 59, 740(1999). 

[21] M. Hasegawa and K. Kaneko, Phys. Rev. C 59, 1449(1999). 

[22] A. Bohr, B. R. Mottelson, and D. Pines, Phys. Rev. 110, 936(1958). 

[23] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 108, 1175(1957). 

[24] M. Manninen, J. Mansikka-aho, H. Nishioka, and Y. Takahashi, Z. Phys. D 31, 
259(1994). 

[25] H. Hakkinen, J. Kolehmainen, M. Koskinen, P. O. Lipas, and M. Manninen, Phys. Rev. 
Lett. 78, 1034(1997). 

[26] W. Satula, J. Dobaczewski, and W. Nazarewicz, Phys. Rev. Lett. 81, 3599(1998). 

[27] R.B. Firestone and V. S. Shirley Table of Isotopes, 8th ed. (Wiley-Interscience New 
York, 1996). 

[28] A. de-Shaht, Phys. Rev. 105, 1528(1957). 

[29] J. Janecke and E. Coniay, Nucl. Phys. A 436, 108(1985). 

[30] E. Comay and J. Janecke, NucL Phys. A 410, 103(1983). 

[31] E. Caurier, P. Zuker, A. Poves and G. Marti'nez-Pinedo, Phys. Rev. C 50, 225(1994). 

[32] J. Duflo and A. P. Zucker, Phys. Rev. C 52, R23(1995). 

[33] J. A. Nolen and J. P. Schiffer, Annu. Rev. NucL Sci. 19, 471(1969). 

[34] M. Hasegawa, Prog. Theor. Phys. Suppl. 132, 177(1998). 

[35] W. Nazarewicz, Nucl. Phys. A 574, 27c(1994). 

[36] W. Nazarewicz et al., Phys. Rev. C 53, 740(1996). 



22 



[37] G. A. Lalazissis and S. Ramman, Phys. Rev C 58, 1467(1998) 



23 



Figure captions 

Fig. 1 Plots of the double differences of binding energies derived from the experimental 

masses in the region A = 16 — 164: (a) 5V^^'{Z, N) as a function o( A = N + Z; (b) 
5V^'^\Z, N) as a function of A = N + Z. The dots stand for even- A nuclei, and the 
crosses for odd-A nuclei. The curve 40 /A is drawn both in (a) and (b). 

Fig. 2 The calculated p — n interaction energies as a function of the valence- neutron number 
nil for (a) the Nb isotopes and (b) the Mo isotopes. The open circles stand for the 
p — n part of the interaction energy E^°'^'^^'^^^, the open squares for the t = p — n 
interaction energy E^^^, and the diamonds for the total p — n energy. 

Fig. 3 The calculated values 5^^^E{Z,N) (open circles) and the experimental values 
6'^^^V{Z, N) (solid squares) as a function of A = N + Z for the Nb, Mo, Tc and 
Pd isotopes. 

Fig. 4 The calculated values 5^'^^E{Z, N) (open circles) and the experimental values 
5'^'^^V{Z, N) (solid squares) as a function of A = N -\- Z for the Mo, Tc, Pd and 
Sn isotopes. 

Fig. 5 The double differences of binding energies for the Ti and Cr isotopes, shown in the 
same manner as Figs. 3 and 4. The solid squares denote the experimental values, 
and the open circles and crosses denote the calculated values with k^ = 1.9 and 
A;0 = 1.9 X {4:8 /A), respectively. 

Fig. 6 The symmetry energy coefficients a{A) in the /7/2 shell region. The calculated and 
experimental values are denoted by the open circles and solid squares, respectively. 
The diamonds represent a{A) = 1.245A in the J = isovector plus J =odd isoscalar 
pairing force model, and the crosses a{A) = 0.295A in the J = isovector pairing 
force model. 
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Fig. 7 Dependence of the experimental odd-even mass differences, As{Z, N) and A^lZ, N) 
on N — Z for nuclei with proton number Z=20-30 and 36. 

Fig. 8 The calculated and experimental odd-even mass differences A^{Z,N) as a function 
oi N ~ Z for the Ca, Ti, Cr and Fe isotopes. 

Fig. 9 Interaction energies of the Pq + QQ + P2 force and VJ^^ in the Cr isotopes. 

Fig. 10 The calculated and experimental odd-even mass differences ^^[Z, N) as a function 
oi N — Z for the Ca, Ti, Cr and Fe isotopes. 

Fig. 11 The two-proton separation energy S2p in the fj/2 shell region. The calculated and 
experimental values are denoted by the open circles and solid squares, respectively. 
The force strength k^ is taken as /c*^ = 1.9 x (48/A). 
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